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Recent X-ray analysis and NMR experiments have shown 
that a'-NaV205 is not a conventional spin-Peierls (SP) com- 
pound because above the transition at 33K it is in a ho- 
mogeneous mixed valent (MV) state and atomic spin chains 
do not exist. Furthermore, thermal expansion experiments 
have identified the existence of a double phase transition. We 
present a theoretical model which qualitatively explains these 
observations. We derive an effective Hamiltonian for the insu- 
lating state which describes the transitions as primary charge 
order leading to linear spin chains in b- direction and a sec- 
ondary SP- transition. We show that this scenario explains 
in a natural way the anomalous BCS- ratio observed in a'- 
NaV 2 5 . 



The discovery of the first inorganic spin-Peierls (SP) 
compound CuGeC>3 JlJ has renewed interest in a phase 
transition originally observed in a few organic compounds 
|||. The spin- chain dimerization in CuGeC>3 below 
Tsp= 14K was investigated in detail using neutron and 
Raman scattering and also by thermal expansion and 
susceptibility measurements. Theoretical understanding 
of the SP- transition as achieved by Pytte In case 
the interchain coupling is sufficiently large a mean field 
HF approximation is possible and then the BCS ratio of 
2A(0)/kT S p=3.52 for the dimerization gap A(0) is ob- 
tained. This is close to the value observed in CuGeC>3 
and in the organic SP- compounds. 

Later on another inorganic 3d- insulator, a'-NaV205 was 
investigated [0. It can be viewed as a layered struc- 
ture with corrugated planes composed of 3d-V bonds ori- 
ented along the segments of zig-zag chains and along the 
'rungs' parallel to the a-axis which connect different zig- 
zag chains (Fig. la). From the temperature dependence 
of the susceptibility the existence of ID spin- chains and 
the opening of an isotropic spin gap below 33K has been 
concluded 0). According to the original X-ray struc- 
ture analysis || of a'-NaV2 05 vanadium atoms exist in 
V 4+ (S=i) as well as V 5+ (S=0) ionization states. The 
V 4+ - species then forms linear AF spin chains along the 
b- axis which has been assumed to be the case at all 
temperatures. Therefore a'-NaV2 05 was interpreted as 
another example of an inorganic SP- system which un- 
dergoes dimerization of the V 4+ - spin chains below 33K. 



Recently it has become clear that this simple view is not 
adequate. Firstly, the BCS- ratio has the anomalously 
high value of 6.44 in a'-NaV205. Furthermore, an im- 
portant observation was made in thermal expansion mea- 
surements H : The phase transition in this compound 
consists of two very close transitions at T cl =33K and 
T C 2=32.7K with first order and continuous character, re- 
spectively. The latter transition is connected with the 
appearance of a large spontaneous strain. This observa- 
tion initiated a new X-ray structure analysis which 
had a very surprising outcome. Contrary to the original 
analysis S a centrosymmetric space group (Pmmn) with 
only one crystallographic vanadium position was found to 
give the best structure refinement at room temperature 
M. This implies that a'-NaV205 is in a homogeneous in- 
sulating MV state with formal valence n=4.5. Therefore 
it cannot be viewed as a system of uniform atomic spin 
chains of V 4+ (S=±) ions above T cl . In addition NMR 
experiments || have confirmed the high temperature ho- 
mogeneous MV state and in addition have proven that 
below the (double) phase transition two incquivalent V- 
sites evolve. 

These findings lead us to the following conjecture: The 
primary order parameter (OP) in a'-NaV205 is a charge 
order of V 4+ /V 5+ below T cl which immediately leads to 
the possibility of a secondary order parameter below T C 2 
connected with the evolution of a spin gap and a large lat- 
tice distortion. Since the latter is characterized by a wave 
vector Q= (j, \, j,) (units of 27r/aetc.) a dimerization 
along the b axis is present and therefore it is possible that 
the secondary OP is mainly of the SP type along the b- 
chain direction. Other possibilities |7],[l(J cannot be ex- 
cluded at present since the precise nature of V 4+ /V 5+ 
charge ordering below T c i is not yet known. However 
inelastic neutron scattering results seem to favor a 
charge order with ID character along b. Irrespective of 
the details a theoretical model for a'-NaV 2 5 must ad- 
dress the following facts: 

(1) at high temperature one has a homogeneous insulat- 
ing MV state. (2) there are two close consecutive phase 
transitions with the primary order parameter driven by 
V 4+ /V 5+ charge ordering. (3) the low- temperature 
phase shows a spin gap with anomalous high BCS- ra- 
tio. (4) there is a lattice distortion including a doubling 
of the unit cell along b- direction. 

In this letter we present and investigate such a model. 
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First we discuss basic features of the electronic structure. 
LDA- calculations Jjjj suggest that one has vanadium d- 
bands with a dispersion of ~ 0.5 eV caused by effective 
hopping (t,t' along the zig-zag chain) and a bonding- 
antibonding gap of 2t ~ leV where t is the effective d-d 
hybridization (Fig. la) via oxygen p-orbitals in the rungs 
between zig-zag chains. Because there is one vanadium- 
3d electron per V-V rung the Fermi level lies in the mid- 
dle of the bonding bands, which means that LDA pre- 
dicts a metallic state. The origin of the insulating state 
of a'-NaV2 05 is therefore not obvious. In Jl^,[l3| it was 
assumed that an on-site Coulomb repulsion U is suffi- 
cient to achieve an insulator. In fact if 2t<U<2£ the 
(half filled) bonding band may acquire a Mott Hubbard 
correlation gap. However in the regime 2t<2i <U< oo 
relevant for a'-NaV205 it is not clear whether an on-site 
U is sufficient to obtain an insulating state or if in addi- 
tion inter-site interactions are needed. 
The essential reason for our introduction of the inter- 
site repulsion V,V,V (Fig.l) however is a different one: 
They are necessary to obtain the observed charge order- 
ing transition irrespective of the nature of the insulat- 
ing ground state. In fact a Hartree Fock calculation 
shows that U> 7,V,V> t,t,V must hold in order 
to achieve complete charge order and moment formation 
at zero temperature. This is also the limit which we con- 
sider here. In this case we take the interaction terms as 
unperturbed Hamiltonian Hq. Its basis states |ni <T ,n2 (T ') 
arc defined by the occupation numbers ni a where 1=1,2 
denotes one of the V- atoms in each rung and er, a' de- 
notes the spin. On the average each rung accomodates 
one d- electron and for an insulating state, irrespective of 
its origin only virtual hopping of electrons between rungs 
is possible. Therefore one can restrict to the one-particle 
(1-p) subspace for each rung which makes it convenient 
to introduce an Ising pseudospin variable r z =n2-ni with 
values ±1 according to whether the 3d electron occupies 
the left (1=1) or right (1=2) V atom of each rung. In 
addition one has the real 3d spin \<5i. In this subspace 
Hq reduces to 
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The kinetic energy Hkin consists of three terms: intra- 
rung hopping (t) and inter-rung hopping (t,f) (Fig. la) 
along the zig-zag chain. Elimination of these processes 
for smalll hopping leads to an effective Hamiltonian ex- 
pressed in terms of r 2 , 07. The former contains "superex- 
change" contributions for the pseudo spins in addition to 
the zeroth- order term of Eq.(l). The latter leads to a 
coupling between pseudo spin (i.e. charge) and spin vari- 
ables. The derivation is indicated for the pseudo- spin 
superexchange originating in t' and V. We divide the 1- 
p sector of each rung with four states into ferromagnetic 
{Po} and antiferromagnetic {Qo} states characterized by 
projectors P =(l+r^rf )/2 and Q =(1-t^t 2 b )/2. The 



energy of {Qo} states is higher by V compared to {Po} 
states. Under the assumption t' < V one obtains an ef- 
fective 2 nd - order contribution in the 1-p subspace given 
by 



H 



(2) 
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It has the same sign (ferromagnetic) as the first- order 
term, i.e., it favors the 3d electron of V 4+ sitting on the 
same side of next neighbor rungs (A,B) (Fig. la). In a 
similar fashion all other kinetic terms (~ t,t') can be 
eliminated. Some of them connect the 1-p sector to the 
2-p sector of doubly occupied rungs. Elimination of such 
terms give either pseudo- spin superexchange along the 
b- axis, or real spin exchange which, however, depends 
on the given pseudo spin or charge configuration. The 
kinetic terms ~ t (intra-rung hopping) do not lead out 
of the 1-p space and therefore may be represented by the 
pseudospin variable rf . One then obtains the effective 
Hamiltonian 



(3) 



(id) 



<l,3> 



(i,j) 



+ J E ^[( 1 + ^)( 1 +^)^ + (l-r:)(l-riK^]. 

By convention denotes nearest- neighbor A,B 

rungs along the [110] direction and <C i,j ~> denotes 
next- nearest neighbor rungs (A, A') or (B,B') which ex- 
ist only along the b- direction. Under the assumption 
that V<U,U the interaction constants obtained from the 
above elimination procedure are given by 
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K' = -V' + — J', K = -V + J, 
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Here the first two terms describe an anisotropic frus- 
trated 2D- Ising model with ferromagnetic n.n and anti- 
ferromagnetic n.n.n interactions in the pseudo- spin vari- 
ables. These terms are responsible for the charge- or- 
dering transition. The third term describes the effect of 
covalency in each rung and it influences the degree as 
well as critical temperature of charge ordering |14f | . In 
the limit t < V' it may also be eliminated and only the 
first two terms of the pseudospin part in Eq.(3) remain, 
but with the replacement K'-» K'+4t 2 /V in Eq.(4). The 
fourth and last terms in Eq.(3) couple charge (r 2 ) and 
spin (^t?) variables in an intricate way: The AF spin ex- 
change ~ J between n.n.n spins along b is only effective 
if electrons are sitting on the same side of the rungs, i.e., 
if the two rungs i,j are charge ordered, which means that 
r*, r| are ferromagnetically aligned. If they occupy oppo- 
site sides, i.e., if r*, r| are antiferromagnetic, there are 
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no virtual hopping processes to doubly occupied states 
possible and hence the spin exchange is zero. When the 
rungs are completely charge ordered along the b-axis the 
last term in Eq.(3) describes an AF S=^ spin chain along 
b. Thus the model of Eq.(3) has the essential ingredients 
to describe the physics of a'-Na 2 V05: 
(1) The high temperature phase corresponds to a homo- 
geneous MV (charge disordered) phase which is insulat- 
ing because inter-rung charge fluctuations are suppressed 
by the U and V interactions. (2) the primary order pa- 
rameter is due to charge order in the vanadium rungs 
which can be described by an anisotropic frustrated 2D 
Ising model for the pseudo spins. (3) Once complete 
charge order has developed, linear S=^ spin chains along 
b (Fig. lb) are formed which leads to the possibility of 
their dimerization as secondary ordering phenomenon. 
In the following an explicit calculation for the order pa- 
rameters is given. As a simple approximation we proceed 
in two steps. The primary OP (r 2 ) is determined by 
the Ising part. In mean- field (MF) approximation the 
charge- ordering temperature would be determined by 
T cl - K(Q) where K(Q)= i(z'K'-zK) (z'=2, z=4) is the 
Fourier transform of the interaction with the charge- or- 
dering wave vector Q c =(0,0) and Q c =(f ,f ) for TM'(+) 
and 'AF'(-) order, respectively. The former corresponds 
to a charge- ordered structure below T c i consisting of 
S=| (V 4+ ) spin chains along b (dotted lines in Fig. lb). 
This is precisely the type of order which was originally 
thought to exist in a'-NaV205 for all temperatures. The 
latter possibility would correspond to a charge ordering 
with V A+ ions filling every second zig-zag chain com- 
pletely. This Ising type model gives a very simplified 
picture of the real charge- order transition in an insula- 
tor where the compensation effect of long range Coulomb 
contributions is important to achieve a low transition 
temperature as in a'-NaV205. In a rudimentary way this 
compensation is included in the Ising terms of Eq.(3) due 
to their opposite signs. The resulting reduction of the or- 
dering temperature can roughly be incorporated by map- 
ping this anisotropic frustrated model to the isotropic 
unfrustrated n.n. Ising model, i.e. by replacing K'— > Ki 
=K'-(z/z')K =K'-K/2 which would be exact in mean- 
field approximation. Note that the direct contribution in 
Ki ~(V'-V) is strongly reduced due to the compensation 
effect. Then the charge- order parameter (r z ) is given by 
Onsager's solution 



(r z ) ={l-sinh(2/3^ 1 )- 4 }i 



(5) 



It is assumed that Ki >0 corresponding to a FM 
Ising OP describing the charge order which leads to the 
straight spin chains along b (Fig. lb). The AF charge or- 
der with isolated zig-zag chains is probably not favored 
JL5| . The charge- order temperature corresponding to 
Eq.(5) is given by T c i=2.27Ki. The order parameter 
is shown in Fig. 2 (dashed line). Its knowledge is essen- 
tial to obtain the effective spin Hamiltonian from Eq.(3). 
thought to be responsible for the secondary ordering at 



T C 2 which opens the spin- excitation gap. Replacing r 2 
— * (t z ) in Eq.(3) we notice the following: The increase in 
(t z ) leads to a rapid quenching of contributions coming 
from the 1=1 chains containing the left V atom of each 
rung (Fig. lb). Therefore only the intra- chain exchange 
contributions from 1=2 chains survive. Consequently the 
inter- chain exchange coupling ~ J' between the 1=1,2 
type chains also becomes ineffective below T c i. Thus the 
Hamiltonian in Eq.(3) describes in a natural way how a 
linear spin- chain system appears below the charge or- 
dering transition. Once it is established it may undergo 
a SP phase transition at T C 2 < T c i leading to a dimer- 
ization order parameter of the linear (1=2) chains along 
b. Its effective ID Hamiltonian is extracted from Eq.(3) 
by defining t z —t z -{t z ) and using 

((1 + ^(1 + 7j)) = [(1 + {r z )f 1] + [1 + (fjf|)]. (6) 

where the first term is the ID- coherent part along b 
and the second term is an incoherent 2D part already 
present above T c i. This separation is equivalent to a 
virtual- crystal approximation for the spin part. Keeping 
only the first coherent part the effective AF spin- chain 
Hamiltonian (Si — \<Ji) then reads 



Hlf f = JT+(T) SiSj, 
{id} 

r + (T) = i[(l + (r,» 2 -l]. 



(7) 



Thus the charge- ordering mechanism leads to a ID 
HAF with an effective temperature dependent n.n. ex- 
change J along the b- chain (this corresponds to n.n.n 
exchange in the plane). The AF- spin chain may now 
exhibit a SP- transition if the coupling to the lattice 
dimerization mode, characterised by a coupling constant 
A, is strong enough. The temperature dependence of the 
dimerization spin gap Asp(T) is obtained from the gap 
equation j|) 



ri(T) 



A| P (T)]i 



x tanh[ 



(r 2 + (T> 



A 2 

^SP 



(T))i 



2kT 



dr. 



(8) 



The cutoff is given by €q— pJ with p= 1.63 0. Note 
that in addition to the intrinsic BCS temperature depen- 
dence there is an additional one caused by the explicit 
appearance of r + (T) in the effective exchange which has 
the following consequences: 

(1) because r + (T)=0 for T>T ci one always has T c2 <T ci 
for the SP- transition. 

(2) the temperature variation of A$p(T) is not BCS- like 
due to the influence of T + (T) and as a consequence the 
BCS- ratio 2Asp(0)/T C 2 will be anomalous. 

We may however define a fictitious SP- transition tem- 
perature T c o through Eq.(8) which would exist if one had 
spin chains with a constant coupling strength jr + (T) = 
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jr + (0) for all temperatures. It is more convenient to use 
T c0 as parameter instead of A. If T c0 > T c i the BCS- 
ratio for the SP transition at T C 2 will become anoma- 
lous. As T c o ^> T c i the physical SP transition at T C 2 
will be very close below the charge- order transition, i.e. 
(T c i- T C 2)/T c i <cl. These findings are clearly illustrated 
in Fig. 2 showing the temperature dependence of both, 
charge ((t 2 ), dashed line) and dimerization (A$p, full 
line) order parameters. We identify T c i with the charge 
order temperature at 33K. Then the only parmeter is 
T c0 /T cl which is fixed to obtain the anomalous BCS- ra- 
tio of 6.44. From this we obtain a SP T c2 of 30.4K or 
(T c i- T c2 )/T c i= 0.08. This is still too high compared to 
the extremely small experimental value of 0.009. How- 
ever this value is determined by the details of the tem- 
perature dependence of (t z ) below T c i which is approxi- 
mated by the 2D Ising model OP. In reality the charge- 
order transition has first- order character Q and this may 
push T C 2 even closer to T c i. 

In conclusion, the present model qualitatively can ac- 
count for the crucial observations in a'-NaV2 05 outlined 
in the beginning. It has the advantage that it gives a nat- 
ural explanation for two seemingly unrelated facts: The 
double phase transition with T c i and T C 2 values close 
to each other and the very anomalous BCS- ratio. As a 
further illustration of this fact we show the dependence 
of this ratio as a function of the T c i- T C 2 - splitting ob- 
tained from our model where T c i/T c o is the parameter 
varied (insert of Fig. 2). Independent of the details of 
the model it has become clear that charge ordering and 
spin (Peierls) gap formation are intricately linked in ol- 
Na2V05 and an important task for future investigations 
is therefore a precise determination of the type of charge 
order. 
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FIG. 1. Top: Charge disordered (mixed valent) structure of 
a'-NaV205 above T c i in the ab- plane. Each 'rung' (dashed 
arrows A,B,A' etc.) contains two vanadium atoms (1,2), one 
in the V 4+ (S=|) state (full circles), the other in the spinless 
V 5+ state (open circles) . The direction of the arrow indicates 
the sign of the pseudo spin (r z ). The configuration shown 
belongs to the 1-p sector for each rung. Hopping matrix el- 
ements t,t',t and inter-site Coulomb interactions V,V',V are 
indicated. They correspond to 2 nd (3.44 A), I s * (3.05 A) and 
3 rd (3.61 A) nearest neighbors. The position of O- and Na- 
atoms is not shown. Bottom: charge ordered (FM-) structure 
of t z below T c i. V 4+ - ions align in linear (2- or 1- type) S=^ 
spin (<?)- chains along b with intra- chain spin exchange J. 
Inter- chain exchange J' is not active in this charge- ordered 
structure. 
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FIG. 2. Temperature dependence of the normalized charge 
({t z ), dashed line) and dimerization (Asp, full line) order 
parameters. Critical temperatures are T c i= 33K and T C 2 = 
30. 4K (arrow) respectively. A coupling parameter A cor- 
responding to T c o= 57K has been used which leads to a 
BCS-ratio 2A 5 p(0)/T c2 =6.47. The insert shows the vari- 
ation of the BCS- ratio with the normalized difference of crit- 
ical temperatures. 
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